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We  seek  to  gain  a  quantitative  understanding  of  the  physical,  chemical,  biological,  and  optical 
dynamics  that  staicture  the  planktonic  ecosystem  on  the  scales  of  the  organisms.  By  quantifying  the 
variability  of  the  organisms  and  their  environment  at  the  scales  that  the  organisms  live  their  lives,  we 
will  gain  insight  into  the  dynamics  that  create  structure  and  influence  the  fiinction  of  the  planktonic 
ecosystem. 


OBJECTIVES 


Our  objectives  in  this  phase  of  the  proposed  work  were  to  develop  image-processing  algorithms  to 
allow  enhanced  and  automated  analysis  of  the  images  obtained  with  our  free-falling  planar  laser 
imaging  fluorometer  (PLIF)  system  (FIDO-<|)).  We  also  interfaced  the  imaging  system  with  an  optical 
nitrate  sensor  and  microstructure  profiler. 

APPROACH 


FIDO-<{>  is  a  free-falling  vehicle  carrying  a  laser,  sensitive  CCD  camera,  computer,  depth/tilt/roll 
sensors,  and  a  CTD/fluorometer/transmissometer  package  (Fig.  1).  The  laser  is  formed  into  a  6.5  mm 
thick  sheet  that  extends  at  a  45°  angle  below  the  vehicle.  Chlorophyll  a  fluorescence  stimulated  by  the 
laser  is  imaged  by  the  camera,  oriented  at  90°  to  the  laser  sheet.  The  imaging  area  is  32x32  cm,  with 
312x312  pm  resolution  per  pixel.  Previous  work  has  shown  that  the  system  can  detect  individual  cells 
at  least  as  small  as  5  pm;  the  fluorescent  organisms  appear  as  individual  bright  points  in  the  images. 

Recently  we  have  modified  the  single-wavelength  PLIF  to  obtain  images  at  2  wavelengths:  685  nm 
(chlorophyll  a  fluorescence),  and  532  nm  (scattered  light  from  the  laser).  We  use  a  filter  wheel  to 
rotate  the  appropriate  filter  in  front  of  the  camera.  Successive  images  of  fluorescence  and  scattering  are 


made  close  enough  in  time  that  the  same  particles  can  be  identified  in  both  images,  and  also  in 
subsequent  pairs  of  images  (which  allows  particle  tracking). 

The  vehicle  falls  at  a  prescribed  rate  (usually  3-10  cm/s),  changing  its  buoyancy  via  pistons  at  a 
predetermined  depth,  and  returning  to  the  surface.  Fluorescence  images  are  acquired  every  2  seconds, 
resulting  in  150-800  images  per  profile,  depending  on  drop  speed.  Data  are  dumped  via  WiFi  while  the 
vehicle  is  at  the  surface,  and  the  vehicle  can  be  re-programmed  (shutter  speed,  drop  speed,  etc.)  for 
subsequent  drops.  Depending  on  battery  power  and  depth,  3-4  profiles  can  be  made  for  each 
deployment,  lasting  about  2  hours. 
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Figure  1.  Left  panel:Arrangement  of  the  laser  sheet,  camera,  and  imaging  plane  in  our  PLIF 
system.  Right  panel:  The  FIDO-<f>  as  it  might  appear  beneath  the  surface,  showing  the  laser  sheet 
extending  at  45°  below  the  vehicle  as  it  falls.  The  FIDO-<p  is  about  3  m  tall,  2  m  wide,  and  weighs 

about  1000  kg  in  air. 


WORK  COMPLETED 

We  deployed  the  single-wavelength  F1D0-<|>  system  at  stations  about  10  km  off  San  Diego,  CA.  The 
bi-spectral  imaging  system  was  deployed  at  stations  in  the  Santa  Barbara  Channel  during  a 
phytoplankton  bloom  that  reached  >50  pg/1  of  chlorophyll  a.  Auxiliary  instruments  on  the  FIDO-  <j) 
frame  included  an  SBE-19  and  SBE-49  CTD,  Wetstar  fluorometer,  ISUS  optical  nitrate  sensor, 
SCAMP  microstructure  profiler,  and  an  RDI  downward- looking  1200  kHz  ADCP.  During  the  Santa 
Barbara  Channel  cruise,  Prof.  Hidekatsu  Yamazaki  (Tokyo  University  of  Marine  Science  and 
technology)  and  his  team  joined  us  and  deployed  the  TurboMAP  turbulence  profiler  while  the  FIDO-<j) 
was  in  the  water.  Additional  data  are  obtained  from  the  shipboard  ADCP  and  CTD/rosette,  and  water 
samples  are  taken  for  phytoplankton  cell  counts  and  size-fractionated  chlorophyll  extractions. 


RESULTS 


The  images  of  chlorophyll  a  fluorescence  clearly  show  individual  phytoplankton  cells,  aggregates,  and 
probably  herbivorous  zooplankton  guts  with  undigested  food.  These  images  provide  a  unique  form  of 
data:  the  two  dimensional  distribution  patterns  of  individual  phytoplankton,  undisturbed,  in  the  ocean. 
From  the  images  we  have  begun  quantifying  spatial  patterns  of  individual  cells,  the  vertical 
distributions  of  the  fluorescent  particle  size-abundance  spectra,  and  the  vertical  variations  of  cell  types. 
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Figure  2.  Particle  area-abundance  profile  from  Drop  14  of  the  FIDO-ty.  Left  panel:  Area- 
abundance  spectra  of  fluorescent  particles.  Black  line  on  all  plots  is  the  total  chlorophyll  a  profile. 
Middle  panel:  Area-abundance  spectra  predicted  from  the  empirical  relationship  of  the  slope  and 
intercept  of  the  spectra  to  the  total  chlorophyll.  Right  panel:  ratio  of  the  observed  and  predicted 
spectra.  Note  the  unusual  layer  of  large  particles  at  35  m.  Horizontal  blue  stripes  show  regions 
where  the  images  were  smeared  due  to  vertical  shear,  and  removed  from  analysis. 


Profiles  of  the  size-abundance  spectra  of  fluorescent  particles  show  a  strong  relation  of  the  abundance 
of  particles  with  the  total  chlorophyll  (Fig.  2).  In  regions  of  higher  total  biomass,  the  size-abundance 
spectrum  had  a  higher  intercept  and  a  steeper  slope.  This  implies  that  there  are  more  of  all  cell  sizes  in 
regions  of  high  biomass;  however  the  largest  cells  are  only  found  in  these  regions.  These  large 
particles  often  formed  thin  (2-3  m  thick)  layers.  We  often  observed  marked  changes  in  the  size- 
abundance  spectral  slope  with  depth.  These  changes  (and  the  implied  changes  in  phytoplankton 
community  structure)  were  often  larger  over  1  m  vertically  than  over  1000  km  horizontally  at  a  single 
depth.  While  the  bulk  trends  of  the  size-abundance  spectra  were  well  predicted  by  empirical 
relationships  with  the  total  chlorophyll,  there  were  often  unusual  layers  that  could  not  have  been 


predicted  from  the  fluorescence  profile  (e.g..  Fig.  2,  35  m).  These  thin  layers  were  not  layers  of  high 
biomass,  but  were  clearly  composed  of  distinct  cell  types,  based  on  the  fluorescent  object  shapes  in  the 
images. 

We  frequently  observed  unpredictable  layers  of  cell  types  (e.g.,  large  cells,  long  cells,  thin  cells,  etc.). 
These  layers  were  quite  apparent  in  the  images,  but  were  invisible  to  a  “regular”  fluorometer  that  just 
measured  total  chlorophyll  a  fluorescence.  There  were  no  obvious  variations  in  total  chlorophyll  that 
would  have  indicated  the  presence  of  these  distinct  layers.  The  ubiquity  of  these  cryptic  layers  shows 
that  the  phytoplankton  community  is  highly  structured,  though  the  bulk  chlorophyll  may  vary 
smoothly.  One  likely  mechanism  for  the  formation  of  these  layers  is  the  vertical  shearing  of  existing 
patches  of  phytoplankton.  The  biomass  (total  chlorophyll)  is  likely  controlled  by  large-scale  external 
forcings  (light  and  nutrients),  while  the  cell  types  making  up  this  biomass  are  a  result  of  stochastic 
factors  such  as  initial  conditions,  grazing,  turbulence,  competition,  etc.  As  horizontally  distinct 
communities  with  the  same  biomass  are  acted  upon  by  the  vertical  shear,  they  become  vertically 
interleaved,  forming  thin,  layered  structures  in  vertical  profiles.  The  processes  generating  horizontal 
patchiness  and  vertical  shear  are  occurring  constantly;  we  predict  that  cryptic  layers  are  a  common  and 
expected  feature  of  marine  planktonic  ecosystems. 

Our  images  give  the  two-dimensional  locations  of  all  the  fluorescent  particles  in  that  volume  of  water. 
We  were  curious  how  the  cells  were  arranged  relative  to  one  another.  Based  on  intuition  and  previous 
work,  we  expected  to  find  cells  clumped  at  some  spatial  scale.  To  quantify  the  distribution  patterns  of 
individual  particles  we  used  the  pair  correlation  statistic.  This  statistic  quantifies  departures  from 
randomness  at  a  range  of  scales.  The  smallest  scale  is  about  twice  the  resolution  of  the  camera  (-0.65 
mm  in  this  case),  while  the  largest  scale  is  about  the  size  of  the  image  (-10  cm).  We  were  surprised  to 
find  that  the  pair-correlation  statistic  showed  that  throughout  water  column,  the  fluorescent  particles  in 
our  images  were  randomly  distributed  on  scales  of  10  cm  and  less.  Spatial  gradients  in  particle 
numbers  are  apparent  over  larger  vertical  scales,  so  there  is  a  transition  from  randomness  to  significant 
gradients  at  about  10  cm.  This  transition  represents  a  changing  balance  of  biological  stratification  vs. 
turbulent  homogenization.  With  no  turbulent  mixing,  the  phytoplankton  would  be  arranged  in  thin, 
stable  layers  according  to  the  light  gradient  and  the  diffusion  of  nutrients.  Turbulence  tends  to  erode 
this  stratification,  leading  to  a  random  distribution  of  cells  at  the  smallest  scales.  We  are  currently 
developing  theory  to  quantify  the  dynamics  linking  biological  stratification  and  turbulence  to  the 
observed  distribution  patterns.  It  is  likely  that  the  vertical  patch  scales  will  be  predictable  based  on 
measurements  of  dissipation  and  density  stratification.  From  this  we  should  be  able  to  predict,  for 
example,  how  far  a  zooplankter  would  have  to  forage  to  find  a  significant  change  in  phytoplankton 
concentration,  given  a  certain  stratification  and  dissipation. 

Mounting  the  ISUS  nitrate  profiler  and  the  SCAMP  microstructure  profiler  on  the  FIDO  has  allowed 
us  to  make  the  first  highly  resolved  estimate  of  vertical  nitrate  flux  (Fig.  3).  While  we  need  to  be 
cautious  in  interpreting  such  limited  data,  it  is  intriguing  that  the  peak  nitrate  flux  occurs  in  a  thin  layer 
at  the  center  of  the  thicker  subsurface  chlorophyll  maximum.  This  is  consistent  with  phytoplankton 
incorporating  the  nitrate,  and  then  becoming  vertically  mixed,  forming  a  thicker  layer  than  the  original 
layer  of  high  nitrate  flux.  We  have  a  great  deal  more  of  this  type  of  data  from  our  recent  cruise  to  the 
Santa  Barbara  Channel,  along  with  better  measurements  of  the  turbulent  dissipation.  These  data  should 
help  us  formulate  and  test  models  of  the  dynamics  supporting  the  deep  chlorophyll  maximum. 
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Figure  3.  Left  panel:  Vertical  profiles  of  the  dissipation  rate  of  turbulent  kinetic  energy  (e),  and 
nitrate  concentration.  Right  panel:  Vertical  profiles  of  chlorophyll  a  concentration  and  the  vertical 
nitrate  flux ,  calculated  as  xf)NOi  idz,  where  the  vertical  diffusivity  k  v  =  0.2eN~2.  Note  that  the  peak 
upward  nitrate  flux  occurs  in  the  subsurface  chlorophyll  maximum. 

IMPACT/APPLICATIONS 

The  results  of  our  research  will  have  a  significant  impact  on  our  understanding  of  the  factors 
controlling  zooplankton  feeding  and  foraging  in  the  ocean,  and  the  role  of  physical  dynamics  in 
structuring  the  planktonic  community  on  scales  of  centimeters  to  100  meters. 

RELATED  PROJECTS 

In  our  NSF-OT1C  sponsored  project  with  Stephen  Monismith  (Stanford  University)  we  have  built  and 
deployed  a  free-falling  stereo  particle  image  velocimetry  system.  By  taking  a  rapid  series  of  images  of 


scattered  light,  followed  by  images  of  fluoresced  light  (both  chlorophyll  a  and  phycoerythrin 
fluorescence),  we  will  are  constructing  three-dimensional  velocity  vectors  with  mm-cm  spatial 
resolution  over  the  imaging  plane. 
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